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Abstract 

Let the reciprocal Newton 'constant' be an apparently non-dynamical Brans-Dicke scalar field 
damped oscillating towards its General Relativistic VEV. We show, without introducing additional 
matter fields or dust, that the corresponding cosmological evolution averagely resembles, in the 
Jordan frame, the familiar dark radiation dark matter dark energy domination sequence. 
The fingerprints of our theory are fine ripples, hopefully testable, in the FRW scale factor; they 
die away at the General Relativity limit. The possibility that the Brans-Dicke scalar also serves as 
the infiaton is favorably examined. 
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ACDM is perhaps the best cosmological model available. Only two free parameters, 
controlling the density of dark matter ^ and dark energy (in the form of a cosmological 
constant A) suffice to parametrize the apparent clash between general relativity (GR) and 
contemporary cosmology. The notion of dark matter is commonly used, but so far we do 
not have any clue whether dark matter is a fact (in which case identifying its constituents is 
a major particle physics challenge) or else an artifact (effectively induced when casting the 
field equations of some beyond GR theory into the traditional Einstein equations format). 
In this Letter we show, without introducing additional matter fields or dust, that the general 
relativistic cosmological interpretation of a damped oscillating^ Newton constant, sliding 
toward its constant VEV, is of a dark radiation dark matter dark energy dominated 
Universe. The cosmological fingerprints of such a theory are fine (suppressed at the GR 
limit) ripples, hopefully testable, in the FRW scale factor. 

The general idea that Nature's fundamental constants may vary in space and time is not 
new. The theoretical motivation for a varying Newton constant, usually accompanied by a 
varying fine structure constant, ranges 

and dilaton/conformalj^ gravity all the way to brane modelsjjand even super-string/M- 



J, and in particular. Zee's 
llH side, however, there is 



theoryjsj cosmology. One may add to this list k-essence models 
broken-symmetric model of gravity. From the experimental 
so far no indication to support a varying Newton constant, with the major cosmological 
constraints coming from CMb[3] and Nucleosynthesis^^. 

Let our starting point be the Brans-Dicke action 

'5 = -T^ / + ^9''''d,(j)d,(j) + V(0) + Crn^ d^x , (1) 
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with (f){x)~^ locally representing the varying Newton 'constant'. The main role of the scalar 
potential V{(f)) will be to allow, by virtue of Zee mechanism [lo|. for the spontaneous emer- 
gence of General Relativity as the BD scalar approaches its constant VEV. In the Jordan 
frame used, the matter Lagrangian Cm is minimally coupled to gravity, and does not couple 
to (pix). In turn, all matter field equations exhibit the standard general relativistic structure, 
and test particles fall freely along geodesies of the (thereby physical) metric g^u- Note that 
a scale transformation g^i, (l)^^gfiu, which is definitely not a general coordinate transfor- 
mation, would give rise to a proper-sign scalar kinetic term plus a non-trivial scalar-matter 
coupling in the so-called Einstein frame. To focus on the physics beyond-GR, however, let 



us momentarily switch off the matter Lagrangian Lm- We will eventually have to switch it 
on, e.g. by adding a tiny amount of ordinary dust, as otherwise it will be impossible to tell 
which metric is actually used by physical rods and clocks. 

On some effective potential grounds (to be clarified soon), we choose the special case 
^BD = 0. Counter intuitively, this will not make (f){x) non-dynamical, and furthermore, 
would not contradict the well-tested GR in the solar system provided jlj] the potential is 
sufficiently confining to render the BD scalar essentially constant. Although the conven- 
tional scalar kinetic term g^'^(t)-^(t)-y has been unconventionally waived away, second order 
derivatives of will still enter the game, owing to the fact that ^—g(l)g^^5TZ^i, is no longer 
a total divergence. Associated with the above action principle, with ubd = ^ enforced and 
in the momentarily absence of Cm-, are the simplified field equations 

By substituting TZ into the trace of the gravitational field equation, the familiar Klein-Gordon 
equation, governed by an effective potential Ve//(0), makes its appearance 

The emergence of such an effective potential actually characterizes the ujbd = case. The 
most general quadratic potential, for example, 

V{(P) = a + 2(3<p + -i(t? , (5) 

gets translated into 

Veff{(l)) = --(2a0 + + const . (6) 
3 

We shall see that, for a suitable range of parameters, such a potential can deliver stable 
attractive gravity, meaning 

(0) ^G-'>0. (7) 

Within the framework of FRW cosmology, with its spatially flat A; = variant favored on 
simplicity (as well as experimental) grounds, the two independent field equations take the 
form 

^H + ^H^ = lvi<P), H + 2H^ = l^-^, (8) 



with H = a/a serving as the Hubble constant. The gravitational field equations can be 
recasted into the Einstein format, that is TZ^^ — \g^uTZ = T^^, thereby defining an effective 
(from GR point of view) conserved energy /momentum tensor. The associated effective 
energy density and pressure are given respectively by 

The GR limit of our generalized Brans-Dicke theory, if exists, can be straight forwardly 
derived, and for the quadratic potential eq.©, one finds 

Clearly, to assure the attractive nature of gravity, one must require a/ [3 < 0. This, together 
with the fact that a non-negative implies {3 > a'')/ (3, furthermore tells us that 7 > (3'^ /a. 
We can finally close the allowed range of parameters, namely 

a>0, /3<0, 7> — >0, (11) 

a 

by insisting on having Veff{(f)) bounded from below. It is important to notice, however, owing 
to the fact that the effective potential is 7- independent, that the background cosmo logical 
constant A > can be arbitrarily fine tuned without affecting Veff{4>). 

On pedagogical grounds, we first discuss the case of a fiat Lorentzian background, and 
then generalize the scope to allow for the cosmological deSitter background. The constraint 
A = can be conveniently parameterized by a = Gcu^G"^, (3 = —Guj"^, and 7 = Gu'^G (this 
mass scale uj has nothing to do with the Brans-Dicke parameter ujbd), so that 

«, ,2 1 
V{<P) = -^(1 - G<pr K//(0) = ^V{<P) . (12) 

Substituting one field equation 

Om = 1 + ^ (if + (13) 

into the other, see eqs.®, leaves us with the second order master equation for the Hubble 
constant 

HH + 3H^H - ^H^ + 2uj^H^ = , (14) 

with LJ being the only mass scale involved. The numerical solution of the latter equations is 
plotted in FigHJ 



H{t) 




FIG. 1: The Hubble constant H{t) and the reciprocal Newton constant (j){t) evolution for A = 0. 
Time averaging over the ripples, justified for large ujt, resembles a General Relativistic matter 
dominated Universe (dashed). 

• At early times {ut <^ 1), the solution of eq. ()14|l is characterized by one relevant dimen- 
sionless constant of integration p 

1 4ujp 



Hit) 



(ujt) 



1/2 



2t 5 



(15) 
(16) 



Eq. ()15p tells us that, at early times, the metric associated with a varying Newton constant 
reminds us of a radiation dominated Universe recovering from a singular Big-Bang. Depend- 
ing on the sign of p, weak initial gravity may be attractive [p > 0) or else repulsive {p < 0). 
It is only the special p = case for which the 0-expansion starts with 

4 



G(j){t) ~ -cuV , 
5 



(17) 



describing strong (small (p), necessarily attractive (positive 0), initial gravity. The latter 
case of preference can be singled out on finiteness grounds when noticing that, as t 0, 
the Ricci scalar and the Lagrangian behave like TZ 12uj'^ (rather than ~ pt^^^"^) and 
ipTZ + V ~ Quj'^G^^ (rather than ~ p^t~^), respectively. 

• At later times {ujt ^ 1), asymptotically approaching GR, the solution exhibits a particular 
1/t-behavior modulated by damped oscillations. This can very neatly numerically verified 



by plotting H ^(t) and noticing that all its minima (corresponding to the maxima of H{t)) 
lie on a straight line. Quite serendipitously, however, is the fact that the slope of this hne is 



J. This observation is fully captured by the analytic expansion 



G(f){t) 



4 cos^ ut 



3t 



1 - 



sin 2u!t + I tan ut + p"" 



2ujt 



2 sin 2ujt 



0(t) 



4 cos ujt 



(18) 
(19) 



3ujt <i)(t) 3t 

where p is a constant of integration (presumably related to the previous p). The correspond- 
ing FRW scale factor a{t) can be integrated now to establish our main result 



sin 2u!t + p 
3ujt 



(20) 



The full numerical solution for a{t) is plotted in Fig|21 (focus at the moment on the A = 
branch). As clearly demonstrated in this Figure, the FRW evolution highly resembles a dark 
matter dominated Universe. Moreover, the non-trivial correlation 



l{t) ~ t^/^ ^ 0^^) ^ Q 



-1 



(21) 



leads to the pleasing conclusion that the emerging matter dominated Universe is in fact 
general relativistic. 




FIG. 2: Rippled FRW scale factor: Dark matter ripples are only 1/t-suppressed, whereas dark 
energy ripples are exponentially suppressed. The ripples die away at the a;t ^ oo general relativity 
limit. 



Two practical questions immediately arise: (i) How large is this u frequency? and 
(ii) Can one experimentally probe the rippled structure? Following the Nyquist- Shannon 
theorem |l5|. we need an experimental sampling frequency u^s > 20; in order to faithfully 
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reconstruct the original signal and avoid aliasing. Obviously, if uj~^ is many orders of mag- 
nitude smaller than the age of the Universe, one would not be able to tell the actual FRW 
scale a(t) from its coarse grained average 



a{t) = 



2T 



t+T 



a{t')dt' 



(22) 



for some T ~ a few uj~^. Tiny deviations from the average would presumably be interpreted 
as measurement errors. In some sense, this reminds us of a tracking oscillating energy 
modelj3] where the potential takes the form of a decaying exponential with small perturba- 
tions superimposed. To appreciate the way the GR limit is approached, consider a physical 
quantity which can be calculated by means of Eq.(j2ni), e.g. the luminosity distance dii^z), 
and compare it with the standard (GR) Hubble plot for a matter dominated Universe (see 
FigEl). 




FIG. 3: Hubble Plot: The luminosity distance d-L^z), plotted as a function of the redshift z, is 
confronted with its corresponding matter-dominated GR rival (dashed). 

Attempting to be a bit more realistic, we now switch on the background cosmological 
constant A. This is done by shifting the potential VA=o{(f)), namely 



VAM<P) = VA=oi<P) + 2GA(P^ . 



(23) 



Notice that Ve//(0) is not sensitive to such a transformation. In turn, V{(f)) and Ve//(0) do 

not share anymore the same absolute minimum, and the vacuum of the theory, reflecting 

dV 

the minimization of VefA(p), is associated now with (7^) = ~ ~4A. In some sense, 

this reminds us of Ref. [l7|, where the constant curvature vacuum solution is deSitter rather 
than Minkowski. At any rate, in the presence of A, Eq. p3|) is generalized into 



G0(t) = 1 + 



2(^^ + U) 



^A 



(24) 



whereas eq. (fT^ acquires a non-homogeneous term 

HH + ?,H^H - -H^ + 2lj'^H^ = -u^A . (25) 
2 3 

The associated numerical solution is plotted in FiglH 

H{t} 
i 1 




FIG. 4: The Hubble constant H{t) and the reciprocal Newton constant (j){t) evolution for A > 0. 
ACDM-like evolution (dashed) is recovered by time averaging over the ripples. 

• At early times (cut ^ 1), as expected, the A = solution eqs. ()15|lf)j) remains qualitatively 
the same. The leading A-modified term in H{t) is only 0{ujH^), while (f){t) is only stretched 
by an overall factor Su"^/ (A + 3c<j^). As before, strong attractive initial gravity is encountered 
for the p = choice, with eq. ()17|) staying intact. 

• At later times {ut ^ 1), however, the A > solution is entirely different from the A = 
solution eqs. p8ll9p . To see what are the possibilities, take H(t) = A/3 + 6(t) and linearize 
eq.(j2Sl) to read 

6 + V3A6 + 4uj^6 = . (26) 

So, it is either damped oscillations with a shifted frequency = u"^ ~ ^A, or else (critical 
or sub-critical) exponential decay. Notice that, in contrast with dark matter ripples, which 
are only 1/t-suppressed, dark energy ripples are exponentially exp (— i-\/3At)-suppressed. 
This is expressed in Fig|21 where one can appreciate the smooth dark matter to dark energy 
transition. The cosmological ACDM-like evolution is recovered (see the dashed lines in FigHJ 
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by time averaging over the ripples. The total number of ripples at the matter dominated 
era is estimated by 

iV~of^^) . (27) 



,7rV3A, 

The structure of the scalar potential V{(f)) is quite arbitrary. As long as a positive Newton 
constant VEV is our only physical requirement, a quadratic potential will do. It is then by 
no means trivial that the entire physical dark radiation — > matter energy domination 
sequence is recovered at no extra cost. But suppose we want more. In particular, one 
would like to incorporate inflation into the theory. Traditionally, a scalar field is anyhow 
invoked for such a purpose, so we are naturally led to a unified inflationary model where 
the infiaton and the reciprocal Newton 'constant' are in fact the one and the same scalar 
field. In which case, to end up with (say) a quartic Higgs-like Ve//(0), one better start, as 
dictated by eq.Q, from a cubic V"(0), namely 

V{(j)) = 4Ao</'(l - G(l))^ + 2GA(j)^ . (28) 

The only constraint imposed on the cubic polynomial, namely V^(0) = 0, which can be 
waived away in case of need, is engineered to give rise to the symmetric (under (p —>■ —(f)) 
double-well effective potential 

VefM = ^(1 - G'<py . (29) 

The coefficient Aq is to be identified as the cosmological constant during the inflationary 
era. 

The non-physical (unstable) yet pedagogical solution (j){t) = 0, that is sitting for ever on 
the top of the Veff{4>) hill, is unexpectedly correlated with a soliton-like solution for H(t) 

Hit) = ^t^nh(^2^?j . (30) 

It represents a smooth pre-Big-BangjlQ^ transition (borrowing string theory terminology, 
one may call it a 'graceful exit'^^) from eternal deflation to eternal inflation. The general 
solution, however, introduces a constant of integration x > 0, such that around 0(0) = we 
have 

^ xJ^t . (31) 
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The so-called 'graceful exit' is now not only a deflation to inflation transition, but further- 
more a novel transition from repulsive {(p < 0) to attractive {(p > 0) gravity. The duration 



of the inflationary period is of order C(— y3/Ao logx), and can be arbitrarily flne tuned by 
adjusting x- The full numerical solution, omitting the interesting < section which is a 
bit beyond the scope of this Letter, is depicted in FigO 



FIG. 5: The Hubble constant H{t) and the reciprocal Newton constant (p{t) evolution governed 
by a double- well Veff{(j)). The novel deflation — > inflation transition is demonstrated. The after 
inflation evolution is the same as in Figll] 

One may flnd it conceptually difficult to go beyond GR, and accept the unconventional 
idea that cosmological dark matter is in fact nothing but a signature of a damped oscillating 
Newton constant. After all, it is easier to phenomeno logically parameterize ACDM than 
struggle with the puzzle where does the cosmological dark matter (and/or dark energy) ac- 
tually come from. In this paper, by minimally generalizing GR, Brans-Dicke style, without 
explicitly introducing any additional matter fields or dust, we have derived a cosmological 
evolution which serendipitously resembles the dark radiation dark matter —>■ dark energy 
domination sequence. The main prediction of our theory is that the FRW scale factor asso- 
ciated with the emerging, necessarily general relativistic, dark matter dominated Universe 
(see Figj21) is superimposed with a rippled structure. To be on the safe side, one can always 
seek shelter at the GR-limit — oo where the ripples die away. But on the other hand, 
one can only hope that the novel mass scale uj is not that large, allowing the ripples to be 




H[t) 



m 



G 
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experimentally tested. Our final remark concerns the time independent radially symmetric 
solution which is still under active investigation [2^. At this stage, we can only confirm that 
the Newton limit has been recovered, with the deviation from the force law being e"'^'' 
Yukawa suppressed. 
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